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abstract 

O ' We investigate the possibility that stellar mass black holes, with masses in 

c/3 the range of 3.8Af0 and 6Mq, respectively, could be in fact quark stars in 

the Color-Flavor-Locked (CFL) phase. Depending on the value of the gap pa- 
rameter, rapidly rotating CFL quark stars can achieve much higher masses 

> 

' than standard neutron stars, thus making them possible stellar mass black 



hole candidates. Moreover, quark stars have a very low luminosity and a com- 
pletely absorbing surface - the infalling matter on the surface of the quark 



00 ........ 

■ star is converted into quark matter. A possibility of distinguishing CFL quark 



stars from stellar mass black holes could be through the study of thin ac- 
cretion disks around rapidly rotating quark stars and Kerr type black holes, 
, respectively. Furthermore, we show that the radiation properties of accretion 

■ disks around black holes and CFL quark stars are also very similar. However, 

strange stars exhibit a low luminosity, but high temperature bremsstrahhmg 
spectrum, which, in combination with the emission properties of the accre- 
tion disk, may be the key signature to differentiate massive strange stars from 
black hole. 

Key words: accretion, accretion discs: black hole physics: dense matter: 
relativity: equation of state. 
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1 INTRODUCTION 



Current physics suggests t hat compact ob.iects wit h mass functions larger than 3-4 solar 



masses must be black holes (IRhoades fc Ruffini 



197^ ). From a theoretical point of view, black 



holes have been considered, in the framework of general relativity, for almost one hundred 
years. However, their first observational evidence was not found until quite recently, during 
the past three decades. In the late 1960s X-ray detectors onboard satellites revolutionized 
astronomy with the discovery of a large population of luminous X-ray sources in the Galaxy. 
Optical counterparts associated with high mass X-ray binaries (HMXBs) have also been 
identified, and one of the first to be identified was the supergiant star HD 226868, associated 
with HMXB Cyg X-1. The determination of the mass function equation / {M^), which relates 



the mass of the compact object with t 



lat o 



angle i indicated that > AMq flBoltoJl975l ). In 1975 the satellite Ariel V detected the X 



the companion star M^. and to the inclination 



ray source A-0620-00, which belongs to a class of X-ray transients (XRT). The mass function 



of the compact star is 3.2 ± 0.2M(7i, which slight 
a stable neutron star ( McClintock fc Remillard 



y exc eeds the maximum mass allowed for 



1986|). In 1989, the X-ray satellite Ginga 



discovered the XRT V404 Cyg. The mass function for this system implies the presence of a 
compact object with mass greater than QMq. Since then, many other stellar mass black hole 
candidates have been found, with seven of them having mass functions in excess of bM^. 
Presently, there are around 20 know n stellar mass blac k holes, like, for example, XTE J1650- 
3.8 solar masses (( 



500, with a mass o: 
3-5 solar masses ( 



Orosz et al. 



Gelino fc Harrison 



20041 ). GRO J0422+32 with a mass between 



20031 ). GRO J1655-40 with a mass of around 6 solar 
masses ( lUdalski et al.ll2003[ ) etc. However, it is estimated that in the Milky Way alone there 
should be at least 1000 dormant black hole XRTs, while the total nu mber of stellar mass 
black holes (isolated and in binaries) could be as large as 100 million (lRomanilll998l ). The 
stellar mass black hole holes have been observed in close binary systems, in which transfer 
of matter from a companion star to the black hole occurs. The energy released in the fall 
heats up the matter to temperatures of several hundred million degrees, and it is radiated in 
X-rays. The black hole therefore is observable in X-rays, whereas the companion star can be 
observed with optical telescopes. The energy released for black holes and neutron stars is of 
the same order of magnitude. Therefore, on the one hand, black holes and neutron stars are 
often difficult to distinguish. On the other hand, since many types of compact objects formed 
from so called exotic matter (boson stars, axion stars, condensed stars etc.) have properties 
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similar to the standard black holes (ITorres 



2002 



Yuan et al. 



20041 ). distinguishing between 



all these classes of compact object is a fundamental challenge for contemporary astronomy 
and astrophysics. There are also important observational evidences for massive black holes, 
like the intermediate-mass black holes (which can be found in the center of globular clusters) 



and the supermassive b 



(IHorowitz fc Teukolsky 



ack h oles in the center of the Milky Way and of the active galaxies 



19991). 



The quark structure of the nucleons, suggested by quantum chromodynamics, indicates 



the possibi li ty of a hadron 



flltoh 



1970 



Bodmer 



1971 



quark 



phase transition at high densities and/or temperatures 



Witten 



19841 ). If the hypothesis of the quark matter is true, 



then some of neutron stars coulc 



(Alcock et al. 


1986 


Cheng et al. 


( 


19981 



Haensel et al 



actu ally be quark stars, built entirely of quark matter 



19861 ). For a general review of quark star properties see 



There are several proposed mechanisms for the formation of quark stars. Quark stars 
are expected to form during the collapse of the core of a massive star, after the supernova 
explosion, as a result of a firs t or second order phase transition, resulting in deconfined 



quark matter (jPai et al.l Il995l ). The proto-neutron star core or the neutron star core is 



a favorable environment for t 



(ICheng et al 



1998 



le co nversion of ordinary matter to strange quark matter 



Chan et al 



20091 ). Another possibility is that some neutron stars in low- 



mass X-ray bi naries can accrete sufficient mass to undergo a phase transition to become 
strange stars (jCheng fc Dailll996h . This mechani sm has also been p roposed as a source of 
radiation emission for cosmological 7-ray bursts (jCheng fc Dailll998l ). 

Most of the investigations of quark star properties have been done within the framework 
of the so-called MIT bag model, with the energy density pc^ and pressure p of a quark- 
gluon plasma related by the equation of state (EOS) p = {p — 4i?)c^/3, where B is the 
difference between the energy density of the perturbative and non-perturbative QCD vac- 
uum (the bag constant). For different values of the bag constant and for different equations 
of state of the quark matter a complete description of the basic astrophysical properties 



mass, radius, eccen tricity, Kepler i an frequency etc.) o f both static and rotating guar 



<i stars 



has been obtained (iWitten 



Gondek-Rosinska et al. 



1984 



2000 



Haensel et al 



1986 



Alcock et al 



1986; 



Dey et al 



1998; 



Harko &: Cheng l2002l ). Quark stars can reach much shorter 



periods than neutron stars, of the order of 0.5 ms. Quark stars could have a radius signif- 
icantly less than that of neutron stars (jCheng et al.lll998l ). r-mode instabilities in rapidly 
rotating quark stars lead to specific signatures in the evolution of pulsars with periods be- 
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low 2.5 ms. If quark matter is absolutely stable, some pulsars, in which the conditions for 
a hadron-quark phase transition are realized, could co nsist of quark matter. Some data on 
pulsar properties are consistent with this assumption ( jMadsenl 120001 ). 

The existence of a large variety of color superconductin g phases of quark matter a t high 



densi t ies has been sugge sted and intensively investigated (lAlford et al. 



1999 



Alford 



2001 



2004 



Alford et al. 



20071 ). At ultra-high density, matter is expected to form a degenerate 



Fermi gas of quarks in which the Cooper pairs of quarks with very high binding energy 
condensate near the Fermi surface (color superconductor). Such a state is significantly more 
bound than ordinary quark matter. It is now widely accepted that at extremely high density 
the ground state of quark matter is the superconducting Color-Flavor-Locked (CFL) phase. 



and tha t this phase q : 



matter (lAlford et al. 



mat ter rather than nuclear matter may be the ground state of hadronic 



20071 ). The existence of the CFL phase can enhance the possibility of the 



existence of a pure stable quark star. The maximum mass of the static quark stars in the CFL 
model is given by M^^^ = IMMq (1 + 5) //B^, where Bqo = B MeV f m~^)/60 MeV fm~^ 
and S = 0.15 (A/100 MeV)^ (60 MeV fm^V^)> with A the gap energy Jflorvath fc Lugones 



20041 ). Hence large maximum masses (M > 4:Mq) can be found for standard values of B and 



A ^ 250 MeV. Rotation can also significantly increase the mass of the CFL quark stars. 
For a discu ssion of the maxiraum m ass of both rotating and non-rotating general relativistic 
objects see IStergioulas et al.l (120031 ). For rotating stars the maximum mass can be of the 
order of Mmax ~ G.IMq (2 x 10"^^ g cm~^ / pm)^^"^ , where pm is the matching density. 

Recent observations suggest that around most of the active galactic nuclei (AGN's) or 
black hole candidates there exist gas clouds surrounding the central compact object, and an 
associate d accretion disc, on a v ariety of scales from a tenth of a parsec to a few hundred 
parsecs (jUrry fc Padovanil Il995l ) . These clouds are assumed to form a geometrically and 
optically thick torus (or warped disc), which absorbs most of the ultraviolet radiation and 
the soft X-rays. The gas exists in either the molecular or the atomic phase. Hence, important 
astrophysical information can be obtained from the observation of the motion of the gas 
streams in the gravitational field of compact objects. 

The determination of the accretion rate for an astrophysical object can give a strong 
evidence for the existence of a surface of the object. A model in which Sgr A*, the 3.7x lO^Mg 







super massive black hole candidate at the Ga 
thermally emitting surface was considered in 



actic center, may be a compact object with a 



Broderick &: NarayanI (120061 ). Given the very 



low quiescent luminosity of Sgr A* in the near-infrared, the existence of a hard surface, even 
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in the limit in which the radius approaches the horizon, places a severe constraint on the 
steady mass accretion rate onto the source, M ^ IO'^'^Mq yr~^. This limit is well below the 
minimum accretion rate needed to power the observed submillimeter luminosity of Sgr A*, 
M ^ 10~^^Mq yr. Thus, from the determination of the accretion rate it follows that Sgr A* 
does not have a surface, that is, it must have an event horizon. Therefore the study of the 
accretion processes by compact objects is a powerful indicator of their physical nature. 

The first comprehensi ve theory of accretio r i disk s around compact general relativistic 
objects was constructed in Shakura fc Sunyaev ( IQTsI ). This theory was e xtended to the gen 



Novikov &: Thorne 



eral re lativistic models of the mass accretion onto rotating black holes in 
( I1973I ). These pioneering works developed thin steady-state accretion disks, where the ac- 
creting matter moves in Keplerian orbits. The hydrodynamical equilibrium in the disk is 
maintained by an efficient cooling mechanism via radiation transport. The photon fiux emit- 
ted by the disk surface was studied under the assumption that the disk emits a black body 



radiation. T^ 
analyzed in 



le properties o 



radian t energ y fiux oyer the thin accretion disks were further 
Page fc Thornd (119741 ) and in iThornd (119741 ). where the effects of the photon 
capture by the hole on the spin evolution were presented as well. In these works the efficiency 
with which black holes convert rest mass into outgoing radiation in the accretion process 
was also computed. 

The emissivity properties of the accretion disks have been recently investigated for dif- 
ferent classes of compact qbiect s, such as rotating and non-rotating boson or fermion stars 
( Torres 2002 : Yuan et al. 20041 ) . as well as for the modified f (R) ty pe theories of grav- 



ity 



Pun et a 



( iHarko et al 



2008al), fo r brane world black holes (jPun et al.ll2008bl ). and for wormholes 



2008 



20091 ). The radiation power per unit area, the temperature of the disk 



and the spectrum of the emitted radiation were given, and compared with the case of a 
Schwarzschild or Kerr-Newman black holes of an equal mass. 

Since quar k stars in the CFL phase can have higher maximum masses than ordinary 
neutron stars (IHorvath fc Lugonesll2004l ). which, in the case of rapidly rotating stars, can 
reach values as high as 6 — 7Mq, the possibility that some stellar mass black holes are 
actually quark stars cannot be excluded a priori. The CFL quark stars, being more massive 
than the presumed neutron-star limit of 3 — 4:Mq, can describe compact objects that do 
not fit into the framework of standard neutron star models. It is the purpose of the present 
paper to propose an observational test of this hypothesis, by considering a comparative study 
of the properties of the thin accretion disks around rapidly rotating CFL quark stars and 
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black holes, respectively. In particular, we consider a comparative study of the basic physical 
parameters describing the disks, like the emitted energy flux, the temperature distribution 
on the surface of the disk, as well as the spectrum of the equilibrium radiation. Due to the 
differences in the exterior geometry, the thermodynamic and electromagnetic properties of 
the disks (energy flux, temperature distribution and equilibrium radiation spectrum) are 
different for these two classes of compact objects, thus giving clear observational signatures, 
which may allow to distinguish CFL quark stars from black holes. 

The present paper is organized as follows. The properties of the general relativistic thin 
accretion disks onto compact objects are briefly described in Section II. In Section III we 
present the metric properties for Kerr black holes as well as for the CFL quark stars. In 
Section IV we consider the radiation flux, spectrum and efficiency of thin accretion disks 
onto Kerr black holes and quark stars. We discuss and conclude our results in Section V. 



2 THIN ACCRETION DISKS ONTO GENERAL RELATIVISTIC 
COMPACT OBJECTS 

Accretion discs are flattened astronomical objects made of rapidly rotating gas which slowly 
spirals onto a central gravitating body, with its gravitational energy degraded to heat. A 
fraction of the heat converts into radiation, which partially escapes, and cools down the 
accretion disc. The only information that we have about accretion disk physics comes from 
this radiation, when it reaches radio, optical and X-ray telescopes, allowing astronomers 
to analyze its electromagnetic spectrum, and its time variability. The efficient cooling via 
the radiation over the disk surface prevents the disk from cumulating the heat generated by 
stresses and dynamical friction. In turn, this equilibrium causes the disk to stabilize its thin 
vertical size. The thin disk has an inner edge at the marginally stable orbit of the compact 
object potential, and the accreting plasma has a Keplerian motion in higher orbits. 

For the gener al relativistic case the theo ry of mass accretion around rotating black holes 



was developed by 



Novikov fc Thornd ( 



973r) . They extended the steady-state thin disk mod- 



els introduced by lShakura fc SunyaevI ( llQTsI ) to the case of the curved space-times, by adopt- 



ing the equatorial approximation for the stationary and axisymmetric geometry. The time- 
and space-like Killing vector fields {d/dtY and [d/dcj))^ describe the symmetry properties of 
this type of space-time, where t and r are the Boyer-Lindquist time and radial coordinates. 
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respectively. A thin accretion disk is an accretion disk such that in cyhndrical coordinates 
(r, 0, z) most of the matter hes close to the radial plane. 

For the thin accretion disk its vertical size (defined along the 2;-axis) is negligible, as 
compared to its horizontal extension (defined along the radial direction r), i.e, the disk 
height H, equal to the maximum half thickness of the disk, is always much smaller than the 
characteristic radius R of the disk, H <^ R. The thin disk is in hydrodynamical equilibrium, 
and the pressure gradient and a vertical entropy gradient in the accreting matter are negli- 
gible. In the steady-state accretion disk models, the mass accretion rate Mq is supposed to 
be constant in time, and the physical quantities of the accreting matter are averaged over a 
characteristic time scale, e.g. At, and over the azimuthal angle = 2tt, for a total period of 
the orbits and for the height H. The plasma moves in Keplerian orbits around the compact 
object, with a rotational velocity Q, and the plasma particles have a specific energy E, and 
specific angular momentum L, which depend only on the radii of the orbits. The particles 
are orbiting with the four- velocity in a disk having an averaged surface density S. The 
accreting matter is modeled by an anisotropic fiuid source, where the rest mass density po 
(the specific internal energy is neglected), the energy fiow vector g'^ and the stress tensor 
^fiu measured in the averaged rest-frame. The energy-momentum tensor describing this 
source takes the form 

T^-^ = pquV + 2M(^g^) + t^" , (1) 

where w^g^ = 0, u^f^" = 0. The four-vectors of the energy and of the angular momentum fiux 
are defined by —E^^ = T^jj{d/dtY and = Tf^^{d/d(j)y , respectively. The four dimensional 
conservation laws of the rest mass, of the energy and of the angular momentum of the 
plasma provide the structure equations of the thin disk. By integrating the equation of the 
rest mass conservation, V^(po^^) = 0, it follows that the time averaged accretion rate Mq 
is independent of the disk radius: Mq = —lirr' Eu^ = const, where a dot represents the 



1974 ). The averaged rest 



derivative with respect to the time coordinate (iPage &: Thornd 
mass density is defined by S(r) = ^^^{pQ)dz, where (po) is the rest mass density averaged 
over At and 27r. The conservation law V ^E^ = of the energy can be written in an integral 
form as 

[MqE -2'nrnW^'']^r = ^T^^/^FE , (2) 

where a comma denotes the derivative with respect to the radial coordinate r. Eq. ([2]) shows 
the balance between the energy transported by the rest mass fiow, the dynamical stresses 
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in the disk, and the energy radiated away from the surface of the disk, respectively. The 
torque in Eq. ([2]) is given by = J^fj{t^^)dz, where (t^^) is the — r component of 
the stress tensor, averaged over At and over a 27i angle. The law of the angular momentum 
conservation, V^J^ = 0, states in its integral form the balance of the three forms of the 
angular momentum transport, 

[MqL - 2TTrW/],r = Any/^FL . (3) 

By eliminating W,//' from Eqs. ([2]) and ([3]), and by applying the universal energy-angular 
momentum relation dE = QdJ for circular geodesic orbits in the form E^r = ^L^r, the flux 
of the radiant energy over the disk can be expressed in terms of the specific energy, angular 



momentum and the angular velocity of the black hol e. Then the flux integral leads to t 



ne 



expression of the energy flux F{r), which is given by (jPage fc Thorndll974j ; lThorndll974j ) 



Fir) = -71^7^%^ r - nL)L,rdr , (4) 



4:7T^/^ {E -ni) 



where the no-torque inner boundary conditions were also prescribed f lPage fc Thorndll974l ). 
This means that the torque vanishes at the inner edge of the disk, since the matter at the 
marginally stable orbit r^s falls freely into the black hole, and cannot exert considerable 
torque on the disk. The latter assumption is valid as long as strong magnetic fields do not 
exist in the plunging region, where matter falls into the hole. 

The presence of magnetic fields would give a non zero torque at the inner edge of the disk. 
If the central object has a high magnetic field, the disk will be truncated at a radius larger 
than that of the last stable orbit for an unmagnetized star. The presence of high magnetic 
fields will also have a significant effect on the energy and angular momentum transport 
between the rotating central object and the accretion disk. A central object surrounded by 
a magnetosphere can also exert a torque on the accretion disk, via the closed magnetic flux 
lines, connecting the disk and the rotating center. Then, the latter injects angular momentum 



into the orbiting plasma, which in turn increases the flux radiated by the disk surface (ILi 



20021 ). In magnetized thin disks, the magnetic energy and magnetic stresses give also a 
contribution to the radiated flux. The gravitational potential is determined by the nature 
of the rotating central object, whereas the presence of the magnetosphere involves various 
magnetohydrodynamic (MHD) effects, complicating the picture. Since the magnetic torque 
i s propor tional to the difference in the angular velocities of the disk and of the magnetic field 



(LI 



20021 ). and the contribution of the magnetic energy and stresses to the disk radiation has 
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the same proportionality factor for magnetized thin accretion disks (IKovacs et al.ll2009l ). this 
simphfication is vahd as long as the plasma co-rotates with the magnetosphere. Numerical 
models of MHD accretion flows indeed show that the rotating plasma can be coupled to the 
magnetic field in turbulent magnetized dis ks where the plasma particles cannot slip along 
the field lines (IMcKinney &: NarayanI 120071 ) . Assuming this type of turbulent accretion we 
remain in the framework of the hydrodynamical description. In the present study we neglect 
all the effects related to magnetic fields, by assuming that the magnetic fields in the central 
object as well as in the disk are very low, and we consider only the pure hydrodynamical 
description of the disk accretion, that is, we focus only on the effects of gravity. 

Once the geometry of the space-time is known, we can derive the time averaged radial 
distribution of photon emission for accretion disks around black holes, and determine the 
efficiency of conversion of the rest mass into outgoing radiation. After obtaining the radial 
dependence of the angular velocity Q, of the specific energy E and of the specific angular 
momentum L of the particles moving on circular orbits around the black holes, respectively, 
we can compute the flux integral given by Eq. (jlj). 

Let us consider an arbitrary stationary and axially symmetric geometry, 

ds'^ = gttdt^ + gtcf>dtd(j) + Qrrdr'^ + geedO"^ + g^^d(f , (5) 

where in the equatorial approximation {\6 — n /2\ ^ 1) the metric functions gu, gt<f>, grr, gee 
and g(j)(j, depend only on the radial coordinate r. The geodesic equations take the form 

dt _ Eg^^ + Lgt^ dcp _ Egt^ + Lgu 



dr gf^ - gttg^4> ' dr gf^ - gug^^ ' 
and 



(6) 



dr^ ^ 



respectively, where r is the affine parameter, and the potential term V{r) is defined by 

V^(^) = -1 + -^^^0<^ + '^ELgt^ + Pgtt 
9i<f, ~ 9tt94)<t> 

For circular orbits in the equatorial plane the conditions V{r) = and Vrif^) = 0, 
respectively, must hold. These conditions give the specific energy E, the specific angular 
momentum L and the angular velocity Q of particles moving on circular orbits around 
spinning general relativistic stars as 

E = - ^ 9u + 9^^ _ 

-gtt - 2gt^Vl - 
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9t4> + g4>4>^ 



-gtt - '2gt4>^ - g^^"^ 



n 



-gt^,r + \/{gt4>,ry - gtt,'. 



(10) 



(11) 



dt g<p<j),T 

The marginally stable orbits r^s around the central object are determined by the condition 



. 



(12) 



V^rr{r)L-^ = 0, which gives 

E"^ g<j)<f),rr + '2ELgt(f)^rr + gtt.rr ~ {g't4> ~ 

By inserting Eqs. (l9l)-( fT0l) into Eq. ( fT2l) . and solving the resulting equation for Tms, we 
obtain the marginally stable orbits, once the metric coefficients gu, gt<f> and g^tp are explicitly 
given. 

The accreting matter in the steady-state thin disk model is supposed to be in thermody- 
namical equilibrium. Therefore the radiation emitted by the disk surface can be considered 
as a perfect black body radiation, where the energy flux is given by F{r) = (TT'^(r) (a is the 
Stefan-Bol tzmann cons tant), and the observed luminosity L (z/) has a redshifted black body 
spectrum (lTorresll2002l ): 



L (U) = ATTd^I (u) = A 



2n 



o /■■'/- ylrd(t)dr 
— cos^ / — — -. 13) 

TTC^ Jn Jo exp [niye/ 1 ) — 1 

Here d is the distance to the source, /(z/) is the Planck distribution function, 7 is the 
disk inclination angle, and rj and Vf indicate the position of the inner and outer edge of the 
disk, respectively. We take = r^s and rf 00, since we expect the flux over the disk 
surface vanishes at r 00 for any kind of general relativistic compact object geometry. The 
emitted frequency is given by Ue = i^(l + z), where the redshift factor can be written as 
1 + Qr sin sin 7 



1 + z 



-gtt - '2^gt<t> 



(14) 



where we have neglected the light bending (ILuminet 



1979 



Bhattacharyya et al. 



2OO1I). 



The flux and the emission spectrum of the accretion disks around compact objects satisfy 
some simple scaling relations, with respect to the linear scaling transformation of the radial 
coordinate, given by r ^ f = r/M, where M is the mass of the compact object. Generally, 
the metric tensor coefficients are invariant with respect of this transformation, while the 
speciflc energy, the angular momentum and the angular velocity transform as E —>■ E, 
L ML and Q Q/M, respectively. The flux scales as F{r) F{r)/M'^, giving the 
simple transformation law of the temperature as T(r) T {r) /M. By also rescaling the 
frequency of the emitted radiation as z/ ^ z/ = z//M, the luminosity of the disk is given by 
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L (z/) — > L (i?) /M. On the other hand, the flux is proportional to the accretion rate Mq, and 
therefore an increase in the accretion rate leads to a linear increase in the radiation emission 
flux from the disk. 

The efficiency e with which the central object converts rest mass into outgoing radiation 
is the other important physical parameter characterizing the properties of the accretion 
disks. The efficiency is defined by the ratio of two rates measured at infinity: the rate of the 
radiation of the energy of the photons escaping from the disk surface to infinity, and the 
rate at which mass-energy is transported to the compact object. If all the emitted photons 
can escape to infinity, the efficiency depends only on the specific energy measured at the 
marginally stable orbit rms, 



e^l - E 



(15) 



There are some solutions for the geometry of the neutron and quark stars where r^s is 
at the surface of the star or even takes values less than the surface radius of the central 
object. In this case the inner edge of the disk touches the surface of the star and plasma 
under the effect of any perturbation due to hydro- or magnetohydrodynamic instabilities in 
the disk will leave the disk and hit the surface. In this case the energy E^, transferred to the 
star from the disk is measured at the radius Rg of the star, and the efficiency takes the form 
e— 1 — -Be, where Ee — E _ . For Schwarzschild black holes the efficiency is about 6%, no 
matter if we consider the photon capture by the black hole, or not. Ignoring the capture of 
radiation by the black hole, e is found to be 42% for rapidly rotating black holes, whereas 
the efficiency is 40% with photon capture in the Kerr potential. For neutron and quark stars 
the efficiency is varying in a broader range, especially if we take into account that E^g and 
Eg can have very different values for different neutron and quark star models. 



3 KERR BLACK HOLES AND CFL QUARK STARS 

In the present Section we briefly review the basic properties of the Kerr black holes and of 
the CFL quark stars. 
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3.1 Kerr black holes 

The Kerr metric, describing a rotating black hole, is given, in the Boyer-Lindquist coordinate 
system, by 

ds^ = - 1 rft^ + 2 a sm^ Odtdcj) + -^rfr^ + 

where E = + cos^ 6, = gl^ — guQij)^ = r"^ + — 2mr, and a is the Kerr parameter 
taking values from to m. In the equatorial plane, the metric components reduce to 

/ 2m\ ma 

9tt = - 1 , gt<j> = 2 , 

V r / r 

oof 2m\ 
9rr = 9^ci, = r +a [1 + — j , 

respectively. 

For a Kerr black hole the geodesic equation Eq. ([7]) for r becomes 

with the effective potential given by 

Veffir) = -1 + 

^^2 _^ ^2^ _^ 2ma^r] + AELmar — (r^ — 2mr) 
These relations may be rewritten in the form 

r^(%)'-V(r). (19) 

with V{r) given by V{r) = r^Aj^Ve//(^), respectively. 

3.2 CFL quark stars 

The chemical composition of neutron stars at densities beyond the nuclear saturation remains 
uncertain, with alternatives ranging from purely nu cleonic compos i tion th rough hyperon 



(18) 



or meson condensates to deconfined quark matter (iPage fc Reddyl l2006l ) . IWittenI (Il984j ) 
proposed that strange quark matter (consisting of up u, down d, and strange s quarks) might 
be the absolute ground state of hadronic matter at all pressures. If strange quark matter is the 
energetically preferred state at very high pressures, but ceases to be energetically preferred 
below some non-zero threshold minimum pressure, then one would expect the existence of so- 
called hybrid stars, with strange quark cores, but with exterior regions composed of normal 
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neutron-star matter at pressures below the threshold value. Therefore the presence of a quark 
core in a neutron star cannot be ruled out, and the hypothesis of the existence of hybrid stars 
may be consistent with the high masses and radii indicated by some recent observations. 
Hybrid stars with q uark cores have been intensively investigated in the physical literature 



Glendenning 



19961 ). In hybrid stars strange quark matter must be in equilibrium with 



normal nuclear matter. The mixed phase can occur, for a reasonable confinement parameter, 
near the normal nuclear saturation density, and goes over into pure quark matter at about 
5 times the saturation density. The onset of mixed and quark phases is compatible with 
some observed classes of neutron stars, but it hinders the occurrence of kaon condensation 



(iGlendenning 



19961 ). However, when quark seeds form in the core of a neutron star, they will 
propagate to the entire star and convert it to the new quark phase. A phase transition occurs 
between the hadronic and quark phase when the pressures and chemical potentials in the 
two phases are equal, Ph = Pq and fih = fiq, where Ph, fJ^h, and Pg, Ug are the pressure s and 
chemical potentials in the hadronic and quark phase, respectively (IGlendenning) Il996l ). The 
change from the metastable neutron matter phase to the stable quark phase occurs as the 
result of the thermodynamical fluctuations in a homogeneous medium, formed of neutrons, in 
which small quantities of the quark phase (called bubbles or nuclei) are randomly generated. 
Once the quark matter is formed inside the neutron star, it will propagate throughout the 
entire star. The combustion of the nuclear matter into quark matter can occur either as a 
slow combustion process, or as a detonation. In both cases, the presence of a quark seed in 
a neutron star will convert t he entire star to a quark star, and the conversion time may be 
of the order of milliseconds (jChan et al.ll2009l ). 

In order to obtain a consistent and realistic physical description of the rotating general 
relativistic quark stars, as a first step we have to adopt the equations of state for the quark 
matter. It is generally agreed today that the color-fiavor-locked (CFL) state is likely to be 
the ground s tate of matter, at least for asymptotic densities, and eve n if the quark masses 



are unequal (lAlford et al 



1999 



Rapp et al 



2000 



Alford et al. 



20071 ). Moreover, the equal 



number of flavors is enforced by symmetry, and electrons are absent, since the mixture is 
automatically neutral. By assuming that the mass of the s quark is not large as compared 
to the chemical potential /x, the thermodynamical po tential of the quark matter in CFL phase 
can be approximated as (iLugones fc Horvatbi2002l ) 



CFL 



47]-2 47r2 



12 In (m,/2/i) 
327r2 



(20) 
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Figure 1. Pressure as a function of density (in a logarithmic scale) for EOS CFL for different values of A. 



where A is the gap energy. With the use of this expression the pressure P of the quark 
matter in the CFL phase can be ob tained as an exphcit function of the energy density e in 



the form (ILugones Sz Horvath 



2002) 



1 , , 
P = - (£ - 45) + 



2X2 



where 



6^ = -a + Ja^ + -772 {e - B) 



(21) 



(22) 



and a = — m^/G + 2A'^/3. In the following the value of the gap energy A considered in each 
case will be also mentioned for the CFL equation of state, so that, for example, CFL800 
represents the CFL EOS with A = 800. For the bag constant B we adopt the value AB = 
4.2 X 10^^ g/cm^, while for the mass of the strange quark we take the value rris = 150 MeV. 

For the considered equation of state of quark matter the pressure- density relation is 
presented in Fig. [TJ 

In order to obtain the emissivity properties of the disk, the metric outside the rotating 
general relativistic stars must be determined. In the present study we calculate the equi- 
lib rium configurations of the rot ating quark stars by using the RNS code, a s intro duced 



m 



Stergioulas fc Friedman! (119951 ) . and discussed in detail in IStergioulas et al 



code was used for the study of different models of rotating neutro n stars in 



(2003). This 



Nozawa et al 



( 119981 ) and for the study of the rapidly rotating strange quark stars (IStergioulas et al.lll999l ). 



The software provides the metric potentials for various types of compact rotating general 
relativistic objects, which can be used to obtain the physical properties of the accretion 
disks. As a first step in our study, by using the RNS code we obtain the mass function of 
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rapidly rotating CFL stars, and determine the maximum mass they can achieve, as well as 
the physical conditions under which this maximum mass can exist. 

The mass distribution of the quark stars in the parameter space {pc, rp/r^) is represented 
in Fig. [21 In our analysis we consider four values of the gap parameter A, namely, A = 
100, 300, 500 and 800 MeV, respectively. The central densities corresponding to these stellar 
configurations are in the ranges of 0.52 — 9 x 10^^ g/cm^ for A = 100 MeV, 0.32 — 9 x 10^^ 
g/cm^ for A = 300 MeV, and 0.3 - 9 x 10^^ g/cm^ for A = 500 MeV and 800 MeV, 
respectively. The lower limits of the central densities are determined by the criterion of the 
positivity of the pressure in the star. For a given value of A we analyze only the configurations 
corresponding to the maximum mass. 

The physical properties of the models with maximal total mass are presented for rotating 
CFL stars in Table [T|, and in Table [2] for static CFL stars, respectively. With the increase 
of the gap parameter the masses of the rotating CFL quark stars increase as well, reaching 
a value as high as 5A6Mq for A = 800 MeV. The increase of the gap parameter A also 
produces a considerable increase in the equatorial radii (from 15 km for A = 100 MeV to 23 
km for A = 800 MeV), and of the angular momentum of the rotating CFL quark stars. In 
each case Re > Rsch, where Rsch is the Schwarzschild radius defined as Rsch = 2GM/c^. For 
A = 100 MeV, Rsch = 7.89 km, smaller than Re = 15 km. For A = 800 MeV, Rsch = 16.18 
km, while Re = 23 km. The angular velocity Q decreases from 10^1/s to 8 x lO'^l/s with the 
increasing values of A. There is also a moderate change in the position of the height 
from surface of the last stable co-rotating circular orbit in the equatorial plane. While for 
A = 100 MeV, /i+ is located on the surface, the last stable co-rotating circular orbits are 
located at 2 and 3 km from the surface of the star in the equatorial plane for A = 300 and 
A = 500 MeV and A = 800 MeV, respectively. 

In the case of the static configurations, in Table [2] we present the physical characteristics 
of the CFL stars with the maximum masses corresponding to the given values of A. The 
maximum masses of the CFL quark stars vary from 1.84Mq for A = 100 MeV to 3.92Mq for 
A = 800 MeV, while the stellar radii change from 10 to 16 km. The corresponding central 
densities pc are in the range of 2.4 X 10^15) g/ cm^ for A = 100 MeV and 9 x 10^^ g/cm^ for 
A = 800 MeV. Both h+ and /i_ have negative values, being located below the surface of the 
star. The localization of h+ strongly influences the energy flux from accretion disks around 
compact objects. 

For all values of the gap energy A the shape of the curves M = M {ve) exhibits the 
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standard form specific to quark stars (jCfieng et al.lll998l ). However, the dependence of the 



mass-radius relation on the central density of the stars is different for the different classes 
of models. For the low mass - low radius limit of these curves (the left bottom part of the 
plots), the lowest possible central density (with its value of 5.2 x 10^^ g/cm^), is considerably 
higher for A = 100 MeV, than the lower limits of pc for A = 300 — 800 MeV, located at 
3 — 3.2 X 10^^ g/cm^. The pc-dependence of the branches running from the maximal masses 
down to the lower mass limit for A = 100, or for higher values of A, is also rather dissimilar. 
Considering, e.g., the case Vp/re = 0.5, for A = 100 MeV and for pc = 6 x 10^^ g/cm^, we 
obtain M = 0.70Mq and Re = 10.72 km, respectively, values that are close to the lower 
limit. The corresponding values in the mass-radius plane for A = 300 MeV are M = 4.7OM0 
and Re = 21.49 km, and even higher for higher values of A, which are located close to the 
maximum of these curves on the mass scale. In order to reach the maximal mass for A = 100 
MeV the central density must be set to 2 x 10^^ g/cm^. For high values of A, the low mass 
branch of these curves can therefore be covered by a much more limited range of central 
densities pc, as compared with the case A = 100 MeV. We conclude that the mass-radius 
relation for A ^ 300 MeV is much more sensitive for the change of the central density than 
it is for A = 100 MeV. Of course, in the latter case the mass-radius range spanned by the 
low mass branch is a factor of two smaller in the mass scale than those belonging to A ^ 300 
MeV, but the maximal change in Re is roughly of the order of 10 — 12 km for any values 
of A. Therefore for A ^ 300 MeV the increase in the equatorial radius is relatively small. 
As compared to nonrotating stars, the effect of rotation is to increase the equatorial radius 
of the star and also to increase the mass that can be sustained at a given central energy 
density. 

One needs to consider stability of the models since models which are unstable with respect 
to radial (or quasi-radial) pulsational modes are not generally of astrophysical interest. A 
quark star with low mass is predominantly self-bound, and is stable to these modes. As 
one moves along the mass-radius curves shown in Fig. [2], proceeding away from the lowest 
mass models, the central density progressively increases. In the absence of rotation, the 
models are stable up to the first maximum in the curve, and unstable thereafter. But with 
rotation the situation becomes more complicated. If such a sequence of models is calculate d 



at constant angular momentum J, then a similar criterion applies (IFriedman et al.lll988l ). 



but our sequences are calculated for constant axis-ratio, and so the criterion does not directly 
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A [MeV] 


100 


300 


500 


800 


Pc [10i5g/cm3] 


1.70000 


0.70000 


0.61500 


0.60000 


M[Mq] 


2.67179 


4.87321 


5.32598 


5.46782 


Mo [Mq] 


3.20367 


6.10475 


6.79655 


6.96692 


Re [km] 


14.8727 


21.9955 


22.4372 


23.2988 


n[io3s-i] 


10.7827 


8.08650 


8.06520 


7.89727 


np[io3s~i] 


10.7541 


8.02614 


8.09665 


7.70780 


T/VK[10-i] 


2.08674 


2.34702 


2.22506 


2.36814 


cJ/GM^ 


6.24670 


21.1302 


24.8740 


26.5042 


/[lO-'^gcm^] 


5.09138 


22.9645 


27.1046 


29.4951 


4>2[10''5gcm2] 


0.87825 


3.95150 


4.50791 


5.10189 


h+ [km] 




2.83756 


3.27146 


3.02038 


h-[km] _ _ _ _ 


ajc/n[io-i] 


7.58826 


8.51031 


8.89081 


8.80415 


re [km] 


9.95880 


12.4564 


11.7922 


12.4806 




0.47000 


0.47000 


0.50000 


0.47000 



Table 1. Physical parameters of rotating maximum mass CFL quark stars for different values of A. Here pc is the central 
density, M is the gravitational mass, Mq is the rest mass. Re is the circumferential radius at the equator, Q is the angular 
velocity, Qp is the angular velocity of a particle in circular orbit at the equator, T/W is the rotational-gravitational energy 
ratio, cJ/GMq is the angular momentum I is the moment of inertia, <I>2 is the mass-quadrupole moment, h+ is the height 
from surface of the last stable co-rotating circular orbit in the equatorial plane, h- is the height from surface of the last stable 
counter-rotating circular orbit in the equatorial plane, ujc/^ is the ratio of the central value of the potential lo to £1, re is the 
coordinate equatorial radius and r^/re is the axes ratio (polar to equatorial). 



apply to them. Roughly speaking, it is likely that stability will continue to change roughly 
at the maxima, but analysis of this is outside the scope of the present paper. 

The differences with respect to the standard quark star models are determined by the 
change of the sign of the quantity 2A^ — in the equation of state of the quark matter 
in CFL phase, Eq. ( 12T|) . If < V^A, the positive correction term will increase the pressure 
in the star, by reducing the numerical value of the bag constant, while for > V2A, the 
negative values of the correction term soften the equation of state, and increase the value of 
the effective bag constant. Hence large A values determ ine a smaller effective bag constant. 



which in turn allows larger stellar masses (1WitteiJll984j ) 



4 DISK ACCRETION ONTO KERR BLACK HOLES AND QUARK STARS 

The metric outside the rotating compact general relativistic stars can be described, in quasi- 
isotropic coordinates, as 

ds^ = -e^+Pdt^ + e^" (df^ + f^dd^) + e^-^f2 sin^ d {dcj) - uodtf , (23) 



18 Kovdcs, Cheng and Harko 



A [McV] 


100 


300 


500 


800 




2.4 


1 


0.9 


0.9 


M [Mq] 


1.84841 


3.41121 


3.84898 


3.92723 


Mo [Mq] 


2.24186 


4.39966 


5.08115 


5.2071 


Re [km] 


9.88262 


15.1998 


16.173 


16.2828 


n [lo-'s-i] 














Up [10*s-l] 


1.5859 


1.12729 


1.09043 


1.0899 


T/W 














cJ/GM^ 














/[gcm^] _ _ _ 


'I>2[gcm2] 














/i+ [km] 








-0.329715 


-4.42561 


[km] 








-0.329715 


-0.442561 
















re [km] 


6.88092 


9.51013 


9.67292 


9.63205 


rp/Ve 


1 


1 


1 


1 



Table 2. Physical parameters of maximum mass static CFL quark stars. 
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Figure 2. Gravitational mass versus equatorial circumferential radius of the rapidly rotating and static CFL quark stars for 
different values of the gap energy A: A = 100 MeV (upper left figure), A = 300 MeV (upper right figure), A = 500 MeV (lower 
left figure), and A = 800 MeV (lower right figure). 
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where the metric potentials 7, p, a and the angular velocity of the stellar fluid relative to the 
local inertial frame u are all functions of the quasi-isotropic radial coordinate f and of the 
polar angle 6. The RNS code produces the metric functions in a quasi-spheroidal coordinate 
system, as functions of the parameter s = f/ (f + fe), where is the equatorial radius of 
the star, which we have converted into Schwarzschild-type coordinates r according to the 
equation r = fexp [(7 — p) /2]. To obtain the radius of the innermost sta ble circular orbits 
rm.K we have used a truncated form of the analytical approximation given by 



Shibata Sz Sasaki 



mm . 



6M [1 - 0.54433g - 0.22619g^ + 0.17989g2 - 0.23002g^+ 

0.26296gg2 - 0.29693g^ + O.UbAGq^Qi) , (24) 



where q = J/M"^ and Q2 = —M2/M^, respectively, and where J is the spin angular momen- 
tum, and M2 is the quadrupole moment. 

Once the metric outside the rotating stars is known, from Eqs. (fTOl) - (fTTI) we obtain the 
angular velocity, the specific energy and the specific angular momentum of the orbiting 
plasma particles in the disk. Then we calculate the integral given by Eq. (jlj), measuring 
the photon flux emitted by the disk surface in thermodynamical equilibrium. From the flux 
we obtain the temperature distribution of the disk, as well as the spectra of the emitted 
radiation. For all our calculations we use an accretion rate of M = 1 x IO^^'^Mq/jt, unless 
specified otherwise. In order to compare the thermal and emission properties of the thin 
accretion disks onto black holes and quark stars, we consider several classes of models in 
which both the Kerr black hole and the quark stars have the same mass and angular velocity. 
The physical parameters of the rotating CFL quark stars used in the present study are 
presented in Table [31 The static case, with lower maximal mass than in the rotating case, is 
also plotted for comparison. 

In Figs. [3]we present the radiated flux distribution over the accretion disk around rotating 
black holes and quark stars, for different values of the total mass M and of the spin parameter 
a* = J/M"^ of the central compact objects. The four plots show that the inner edge of the 
accretion disk is located at a lower radius, and the amplitude of the flux radiated from the 
disk surface increases with increasing total mass and spin parameter for both type of objects. 
However, the inner radius is greater, and, consequently, the values of F{r) are smaller 
for quark stars than for black holes in each configuration of M and a*. For comparison 
we also plotted the flux profiles of the accretion disk for non-rotating CFL stars with the 
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A [MeV] 


100 


300 


500 


800 


Pc [IQlSg/cm^] 


1.50000 


1.75000 


0.61500 


0.60000 


M[Mq] 


2.01099 


3.80791 


5.04850 


5.46782 


Mo [M@] 


2.40633 


4.76377 


6.43507 


6.96692 


Re [km] 


11.4128 


16.1807 


21.1469 


23.2988 




7.70653 


7.71710 


7.69576 


7.89727 


np[10*s-i] 


1.33266 


1.06480 


8.61193 


7.70780 


r/w[io-2] 


8.29980 


11.5253 


20.6856 


23.6814 


cJ/GM^ 


2.22209 


9.45230 


21.3040 


26.5042 


/lion's gcm^] 


2.53405 


10.7646 


24.3289 


29.4951 


*2[10"'5gcm2] 


2.13271 


10.2406 


37.3391 


51.0189 


/i+ [km] 


0.00000 






3.02038 


?i_[km] 




-0.155693 


0.00859 


4.78999 


aJc/n[10-i] 


6.44529 


8.76271 


8.61004 


8.80415 


re [km] 


7.99607 


9.19234 


11.1606 


12.4806 


rp/re 


0.78000 


0.73500 


0.57000 


0.47000 



Table 3. Physical parameters of the rotating CFL quark stars used for the calculation of the electromagnetic properties of the 
thin accretion disks. 



maximum equilibrium mass M static- Since the geometry outside the static star is given by 
the Schwarzschild metric, the properties of the flux profiles for the static star coincide with 
those for static black holes. The inner edge of the disk is located at the radius r = 6M, and 
the maximal flux is lower than the one of the rotating cases, in spite of the fact that Mstatic 
is smaller than M. The latter indicates that rotation has a considerable effect in producing 
a higher disk radiation flux, as compared to the disk luminosities for the static cases. All 
these distinctive features for the rotating and static cases can be seen in the profiles of the 
disk temperature as well, which are presented in Figs. |H 

The characteristics of the disk spectra around black holes and quark stars also exhibit 
significant differences. For black holes, the maxima of the spectra are located at higher 
frequencies, and reach higher values than those for quark stars. In case of the static CFL 
quark stars with the maximal total mass, both the spectral amplitude, and the frequency 
at which the maximum is located, have lower values than those for the rotating cases (with 
higher maximal masses). 
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Figure 3. The time-averaged flux radiated by a thin accretion disk around rotating black holes and quark stars with the same 
total mass M and spin parameter a, . The static case with lower maximal mass is also plotted for comparison. 





M [Mq] 


2 


3.8 


5 


5.5 




a. 


0.5494. 


0.6519 


0.8358 


0.8865 


EOS CFL 


fin [km] 

€ 


14.64 
0.0762 


21.97 
0.0921 


24.84 
0.1130 


26.31 
0.1228 


BH 


Tin [km] 
e 


11.99 
0.0863 


20.30 
0.0971 


20.25 
0.1315 


19.46 
0.1495 



Table 4. The radius of the inner disk edge ri„ and the efficiency e for quark stars with EOS CFL and Kerr black holes with 
the same total mass and spin parameter. 

5 DISCUSSIONS AND FINAL REMARKS 

Table m shows the conversion efficiency e of the accreted mass into radiation for both rotating 
black holes and quark stars. For a given configuration of total mass M and spin parameter 
a*, e is somewhat higher in the accretion process driven by black holes, as compared to 
the values obtained for rotating CFL quark stars. This means that rotating black holes can 
always convert more efficiently mass into radiation than CFL type rotating quark stars can 
do. 

As shown by the flux integral in Eq. (j4]), and the explicit expressions of the specific 
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Figure 4. The temperature profiles of the thin accretion disk around rotating black holes and quark stars with the same total 
mass M and spin parameter a* . The static configuration is also presented. 



energy, specific angular momentum and angular velocity given by Eqs. - flTT]) . the rather 
different characteristics of the radial flux distribution over the accretion disk, of the disk 
spectra, and of the conversion efficiency, respectively, are due to the differences between the 
metric potentials of the black holes and of the quark stars, respectively. Even if the total mass 
and the angular velocity are the same for each type of the rotating central object, producing 
similar radial profiles for Q, E and L, the radiation properties of the accretion disks around 
these objects exhibit large differences. The reason is that the proper volume, and in turn 
the function y/—g, used in the calculation of the flux integral, is highly dependent on the 
behavior of the metric component grr = {dr /dr^grr- The latter contains the r-derivatives of 
the metric functions p(r) and 7(r), respectively, via the coordinate transformation between 
r and r, which are extremely sensitive to the slopes of p(r) and of 7(r), respectively. As a 
result, though the inner edges of the disks are located at almost the same radii in the flux and 
temperature plots presented in Section HI the amplitudes of the energy fluxes emerging from 
the disk surfaces around black holes and CFL quark stars are significantly different. Similar 
differences also appear in the disk spectra of quark stars and black holes, respectively. 
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Figure 5. The disk spectra for rotating black holes and quark stars with the same total mass M and spin parameter a*. The 
static case 



The Coulomb barrier at the quark surface of a hot strange quark star may also be a 
powerful source of e~^e~ pairs, which are created in the extremely strong electric field of the 
barrier. At surface temperatures o f around 10 ^^ K, the luminosity of t he outflowing plasma 
m ay be of the order ~ 10^^ ergs~^ ( lUsov 



1998 



by 



Harko fc Cheng]|2006l ). Moreover, as shown 



Page &: Usovl (|2002| ). for about one day for normal quark matter and for up to a hundred 



years for superconducting quark matter, the thermal luminosity from the star surface, due 
to both photon emission and e^e~ pair production may be orders of magnitude higher than 
the Eddington limit. 

Photon e missivity is th e basic parameter for determining macroscopic properties of stellar 

(119861 ) have shown that, because of very high plasma frequency Up 



type objects. 



Alcock et al. 



near the strange quark matter edge, photon emissivity of strange quark matter is very low. 
Propagation of electromagnetic waves of frequencies lower than Up is exponentially damped. 
Therefore, only photons produced just below the surface with momenta pointing outwards 
can leave strange quark matter. For temperatures T « Ep/uj, where Ep ^ 23 MeV is 
the characteristic transverse plasmon cutoff energy, the equilibrium photon emissivity of 
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strange quark matter is negligible small, compared to the black body one. The spectrum 
of equilibrium photons is very hard, with hu > 20 MeV. The problem of the soft photon 
emissivity of quark ni atter at the surface o f stran ge quark stars has been considered in 



Cheng fc Harkd (120031 ) and iHarko fc ChengI (120051 ) . By taking into account the Landau- 



Pomeranchuk-Migdal effect and the absorption of the radiation in the external electron 
layer, the emissivity of the quark matter can be six orders of magnitude lower than the 
equilibrium black body radiation. However, despite these very specific signatures for quark 
stars, a definite method for discriminating them with respect to the neutron stars is still 
missing. Quark stars have a very low photon emissivity, which make them very difficult to 
directly identify them observationally. 

Unpaired strange quark matter in bulk contains u, d, s quarks as well as electrons to 
maintain charge neutrality. The chemical balance is maintained by weak interactions and 
neutrinos assumed to escape from system. If strange quark matter is in the CFL phase, in 



which quarks of all flavors and colors near the Fermi surface fo rm pairs, an 



of flavors is enforced by symmetry, and the mixture is neutral (lAlford et al.l 120071 ). There 



equa 



number 



fore CFL quark stars do not have an el ectrosphere, and the main energy loss m echanism 



is through quark-quark bremsstrahlung (jCheng fc Harko 



2003 



Chan et al.l 



20091). The lu- 



minosity Ls of the accreting quark star is given by Ls = rjMc^, where for the energy 
conversion rate rj we can take the value t] ^ 0.1. For an accretion rate M = 10~^^Mo/year, 
the luminosity of the star is Ls ~ 5.78 x 10'^^ erg/s, while the luminosity Ld of the disk 
for the same accretion rate is, depending on the value of the gap energy A, in the range 
Ld ^ 2 X 10^^ — 4 X 10^^ erg/s. Therefore the luminosities of the disk and of the quark 
star are very close, having approximately the same values. On the other hand, the surface 
temperature of the central accreting quark star is given by T.g = {riMc^ I AnR^ea] ^ , wher e 



Chan et al 



20091 ). 



R is the radius of the star, and e = 10^^ — 10~^ ( ICheng &: Harkd 12003 : 
For e = 10~^, the surface temperature T5 of the star is of the order of T5 3 x 10^ K, 
while the corresponding temperature of the disk is around T^^ ^ 6 x 10^ — 9 x 10^ K. 
Therefore the temperature of the star is much larger than that of the disk, Ts » T^. 
Furthermore, while the spectral distribution {v) of the radiation coming from the disk 
is a standard black body di stribution, the surfac e radiation of the quark star has a typical 
bremsstrahlung spectrum ( Harko fc Cheng) l2005l ). Therefore, the observational detection of 
a bremsstrahlung type electromagnetic emission from the central object, in combination 
with the emission properties of the accretion disk, could provide a distinct signature for the 
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presence of a high mass and extremely low luminosity quark star in a black hole candidate 
system. 

In conclusion, the observational study of the thin accretion disks around rapidly rotating 
compact objects can provide a powerful tool in distinguishing between standard black holes 
and stars with exotic equations of state, that have underwent, for example, a phase transition 
from the neutron phase to a quark phase, as well as for discriminating between the different 

phases of the dense matter. 
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